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Abstract 

The dimerized one-dimensional Hubbard model is studied in the framework of lattice density- 
functional theory (LDFT). The single-particle density matrix with respect to the lattice sites is 
considered as basic variable. The corresponding interaction-energy functional VF[7jj] is defined by 
Levy's constrained search. Exact numerical results are obtained for VF(7i2 , 723) where 712 = 7i,i+i 
for i odd and 723 = 7j,j+i for i even are the nearest-neighbor density-matrix elements along the 
chain. The domain of representability of 7jj and the functional dependence of W^(7i2 ) 723) are 
analyzed. A simple, explicit approximation to VF(7i2,723) is proposed, which is derived from 
scaling properties of W, exact dimer results, and known limits. Using this approximation, LDFT 
is applied to determine ground-state properties and charge-excitation gaps of finite and infinite 
dimerized chains as a function of the Coulomb-repulsion strength U jt and of the alternation bt of 
the hopping integrals tij {tij = t±5t). The accuracy of the method is demonstrated by comparison 
with available exact solutions and accurate numerical calculations. Goals and limitations of the 
present approach are discussed particularly concerning its ability to describe the crossover from 
weak to strong electron correlations. 

PACS numbers: 71.15.Mb, 71.10.Fd 
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I. INTRODUCTION 



Hohenberg and Kohn (HK) replaced the wave function by the electronic density p{r) as 
the fundamental variable of the many-body problem and thereby achieved a crucial break- 
through in the theoretical description of the electronic properties of matter.0 Since then, 
density functional theory (DFT) has been the subject of remarkable developments. Formal 
improvements, extensions, and uncountable successful applications to a large variety of prob- 
lems have made of this theory the most efficient, albeit not infallible, method of determining 
physical and chemical properties of matter from first principles.i'i Density functional (DP) 
calculations are usually based on the Kohn-Sham (KS) scheme which reduces the correlated 
A^-particle problem to the solution of a set of self-consistent single-particle equations.^ While 
this transformation is formally exact, the form of the interaction-energy functional iy[p(r)] 
involved in the KS equations is not known explicitly. Practical implementations of DPT 
always require approximations to iy[p(r)], or equivalently to the exchange and correlation 
(XC) functional Exc[p{r)], on which the quality of the results depends crucially. Therefore, 
understanding the functional dependence of IV[p(r)] and exploring new ways of improving 
its approximations are central to the development of DP methods. 

The most extensively used forms for 14^[p(r)] are presently the local density approxima- 
tion (LDA),@ its spin-polarized version or local spin-density approximation (LSDA)i and 
the gradient corrected extensions,! which were originally derived from exact results for the 
homogeneous electron gas. Despite an unparalleled success in the most diverse areas, the 
LDA-based approach fails systematically in accounting for phenomena that involve strong 
electron-correlation effects as observed, for example, in Mott insulators, heavy-fermion ma- 
terials, or high-Tc superconductors. These systems are usually described in the framework of 
parametrized lattice Hamiltonians such as Anderson,^ Hubbard,! Pariser-Parr-Pople,i and 
related models which focus on the most relevant electron dynamics at low energies. How- 
ever, even with simplified model interactions and a minimal number of orbitals per atom, a 
detailed understanding of the electronic properties in the strongly correlated limit remains 
a serious theoretical challenge. Exact results are rare or numerically very demanding and a 
variety elaborate many-body techniques have been specifically developed in order to study 
this problem.0 Being in principle an exact theory, the limitations of the DP approach have 
to be ascribed to the approximations used for exchange and correlation and not to the un- 
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derlying formalism. It is therefore very interesting to extend the range of apphcabihty of 
DFT to the many-body lattice models that describe the physics of strongly correlated sys- 
tems. Moreover, the development of lattice density-functional theory (LDFT) constitutes 
an intrinsically inhomogeneous approach and provides a true alternative to the LSDA and 
related gradient-corrected methods. Thus, studies on simple models can open new insights 
into the properties of W that should also be useful for future extensions to more realistic 
Hamiltonians and first principles calculations. 

Several physical problems have been already investigated by applying the concepts of 
DFT to lattice models, for example, the band-gap problem in semiconductors,^ the role 
of off-diagonal elements of the density matrix and the non-interacting v representability 
in strongly correlated systems,@ or the development of energy functionals of the density 
matrix with applications to Hubbard and Anderson models.lli In previous work^lEl we have 
considered a density-matrix functional theory of many-body lattice models, that is analogous 
to Gilbert's approach in the continuum,lll and applied it to the Hubbard Hamiltonian with 
uniform nearest-neighbor (NN) hopping integrals tij = t. The interaction energy W of 
the Hubbard model has been calculated exactly as a function of the density matrix for 
various periodic lattices having jij = 712 for all nearest neighbors i and j. An analysis 
of the functional dependence of ^^(712) for different band fillings and lattice structures 
revealed very interesting scaling properties.lii On this basis, a simple general approximation 
to 1^(712) has been derived which yields a remarkable agreement with available exact results 
in ID systems and which predicts successfully the ground-state energy and charge-excitation 
gap of the 2D Hubbard model in the complete range of interaction strength.0 This shows 
that DFT with an appropriate approximation to W is an efficient tool for determining the 
electronic properties of many-body lattice models. 

The purpose of this paper is to extend the method by allowing for alternations of the 
density-matrix elements 7jj between nearest neighbors in order to study the dimerized ID 
Hubbard model. This problem has motivated a considerable research activity in past years, 
particularly concerning the role of electron correlations in the dimerization of polymer chains 
like polyacetylene.0 In this context two qualitatively different regimes may be distinguished 
depending on the relative importance of the intra-atomic Coulomb repulsion U and the NN 
hopping integral t. On the one side, for small U/t, the dimerization can be regarded as a 
bond-order wave that opens a gap at the Fermi surface of the one- dimensional (ID) single- 
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particle band-structure (Peierls distortion). On the other side, for large U /t, local charge 
fluctuations are severely reduced and the low-energy properties are dominated by antiferro- 
magnetic (AF) correlations between spin degrees of freedom. In this case the dimerization 
can be regarded as an alternation of the strength of AF correlations along the chain or 
spin-Peierls state. One of our aims is to analyze the differences between these two types of 
behaviors in the framework of LDFT. The properties of dimerized chains are also very inter- 
esting from a purely methodological point of view. They provide in fact a simple, physically 
motivated means of exploring the functional dependence of Vr[7ij] by including additional 
degrees of freedom, thereby allowing for a larger flexibility. Moreover, several exact results 
are available to compare with (e.g., Bethe-Ansatz solution for the non-dimerized Hubbard 
chain, flnite-ring Lanczos diagonalizations, or density-matrix renormalization-group calcu- 
lations) which allow to quantify the accuracy of the flnal results. 

The remainder of the paper is organized as follows. In Sec. H the main steps in the for- 
mulation of LDFT are briefly recalled. The properties of the interaction-energy functional 
W of the dimerized Hubbard model are discussed in Sec. |T|. The domain of representability 
of and the scaling behavior of W are investigated. A simple explicit approximation to W 
is derived, which is appropriate for direct calculations. Sec. ^ is concerned with applica- 
tions. The ground-state energy and the charge-excitation gap of flnite and inflnite dimerized 
chains are determined as a function of Coulomb-repulsion strength U /t and hopping-integral 
alternation 5t. The LDFT results are contrasted with accurate numerical solutions in order 
to discuss goals and limitations of the present approach. Finally, Sec. summarizes the 
main conclusions and points out some future perspectives. 

II. DENSITY FUNCTIONAL THEORY ON A LATTICE 

We consider the many-body Hamiltonian 

ijcF klmn 

where c|g. (cjo-) is the usual creation (annihilation) operator for an electron with spin a at site 
or orbital %. The hopping integrals tij deflne the lattice (e.g., ID chains, square or triangular 
2D lattices) and the range of single-particle interactions (e.g., up to flrst or second neighbors). 
From the ah initio perspective tjj is given by the external potential Vey±{f) and by the choice 
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of the basis. Vkinm defines the type of many-body interactions which may be repulsive 
(Coulomb like) or attractive (in order to simulate electronic pairing) and which are usually 
approximated as short ranged (e.g., intra-atomic). Eq. (|I]) is mainly used in this section 
to present the general formulation which can then be applied to various specific models 
by simplifying the interactions. A particularly relevant example, to be considered in some 
detail in following sections, is the single-band Hubbard Hamiltonian with NN hoppings.!! 
The non-dimerized form of this model is obtained from Eq. (|l|) by setting tij = —t for i and 
j NN's, tij = otherwise, and Vkinm = U5ki5nmhn- 

The hopping matrix tij plays the role given in conventional DFT to the external po- 
tential Vextij)- Consequently, the single-particle density matrix 7^ between lattice sites 
replaces the continuum density p(r) as basic variable. The situation is similar to the density- 
matrix functional theory proposed by Gilbert for the study of non-local pseudo-potentials 
K!xt(^, ^)-00l^ The ground-state energy Egg and density matrix 'yfj are determined by 
minimizing the energy functional 

EYi,,]=EKb.j\ + WYl,j\ (2) 
with respect to 7ij. E[^ij\ is defined for all density matrices that can be written as 

= E 7ii<x = E 1 1 ^) (3) 

for all i and j, where is an A^-particle state. In other words, 7jj derives from a physical 
state and is said to be pure-state A^-representable.0 The first term in Eq. (^ is given by 

EK = J2^ij7ij. (4) 

It represents the kinetic energy associated with the electronic motion in the lattice and 
includes all single-particle contributions. Notice that Eq. (^ yields the exact kinetic energy 
and that no corrections on Ex have to be included in other parts of the energy functional as 
in the KS approach. The second term in Eq. (^ is the interaction-energy functional given 
byil 
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Vr[7jj] = min 



klmn 



(5) 



where the minimization implies a search over all A^-particles states |^[7ij]) that satisfy 

(^[7^,]| E49. I^bi,]) =7i.- (6) 



for all i and j. W['yij] represents the minimum value of the interaction energy compatible 
with a given density matrix 7,^. It is often expressed in terms of the Hartree-Fock energy 

Eupllij] = o X! ^^i*^' ilijajkla' - Saa'^ilalkja) (7) 

^ ijkl 

and the correlation energy Ed'jij] as 

W[j.,] = EnFhj] + Ech,] . (8) 

W and Eq are universal functionals of jij in the sense that they are independent of tij, i.e., 
of the system under study. They depend on the considered interactions or model, as defined 
by Vkimni on the number of electrons N^., and on the structure of the many-body Hilbert 
space, as given by and the number of orbitals or sites Na- 

E[^] is minimized by expressing = 7jjt+7yi terms of the eigenvalues rjka (occupation 
numbers) and eigenvectors Uika (natural orbitals) as 

lij'y = J2 ^ika Vka U*^^ . (9) 
k 

Lagrange multipliers /i and \ka {^ka = ^ka/Vka) are introduced in order to impose the 
usual constraints J^kaVka = and J2i l^ifco-P = 1- Derivation with respect to m*^^ and rik^. 



(0 < rjka < 1) yields the eigenvalue equationsli^'lla 

%<7 2^ \ + I UiUa = SkaUjka (10) 

with the following conditions relating rjka and Eka- 

Ska < if rika = l , (11) 

Ska = IJ- if < r^fc^ < 1 , (12) 

and 

Eka > /i if r]ka = . (13) 

Self-consistency is implied by the dependence of dW/d'^ij on rjka and Uika- Eqs. (|lO|)-(p!3D 
hold exactly in all interaction regimes. They are analogous to well-known results of density- 
matrix functional theory in the continuum.0 However, notice the difference with the KS-like 
approach considered in Ref.0, which assumes non-interacting f-representability, and where 
only integer occupations are allowed. 



The importance of fractional natural-orbital occupations has already been stressed in 
previous density- matrix functional studies in the continuumS In fact, in the case of models 
one observes that < rjka < 1 for all ka except in very special situations like the uncorrelated 
limit {Vkimn = 0) or the fully-polarized ferromagnetic state in the Hubbard model. This can 
be understood from perturbation-theory arguments — none of the r^^o- should be a good 
quantum number for Vkimn 7^ — and has been explicitly demonstrated in exact solutions 
for finite systems or for the ID Hubbard chain.ii Therefore, the case (|l^ is the only relevant 
one in general. All Eka in Eq. ([To|) must be degenerate or equivalently 

tij + dW/d-fij^ = 5ij II . (14) 

Clearly, approximations of W in terms of the diagonal 7^ alone can never yield such a 
behavior. Given a self-consistent scheme that implements the variational principle, the 
challenge remains to find good approximations to W^[7jj] that are simple enough to be 
applied in practical calculations. 



III. INTERACTION-ENERGY FUNCTIONAL 

In order to determine W^[7jj] from Eq. we look for the extremes of 

klmn 
<t<t' 

+ T.^^3 ((^iE4c,.i^)-7.,) (15) 

hi <^ 

with respect to |\E'). Lagrange multipliers e and \ij have been introduced to enforce the 
normalization of |\E') and the representability of 7jj as required by Eq. (|^). Derivation with 
respect to (^E'l, £ and \ij yields the eigenvalue equations^ 

E 4c,. 1^) + E ^ki^n cU^ycnaCi^ \m) = e\^) , (le) 

klmn 



and the auxiliary conditions (\I'|\I') = 1 and 7jj = (\E'| X^o- 4cjo-|^)- The Lagrange multipliers 
\ij play the role of hopping integrals to be chosen in order that yields the given 7jj. The 
pure-state representability of 7jj ensures that there is always a solution. The subset of 7ij 
that can be represented by a ground-state of Eq. (|TB|) for some Xij is the physically relevant 
one, since it necessarily includes the absolute minimum -^fj of E[^ij\. Nevertheless, it should 



be noted that pure-state representable may be considered that can only be described by 
excited states or by linear combinations of eigenstates of Eq. ([T6|) .B 

The general functional Vr[7ij], valid for all lattice structures and for all types of hy- 
bridizations, can be simplified at the expense of universality if the hopping integrals are 
short ranged. For example, if only NN hoppings are considered, the kinetic energy is 
independent of the density-matrix elements between sites that are not NN's. Therefore, the 
constrained search in Eq. (|^) may be restricted to the |\&[7jj]) that satisfy Eq. (|]) only for 
i = j and for NN ij. This reduces significantly the number of variables in Vr[7ij] and renders 
the determination and interpretation of the functional dependence far simpler. 

In Sec. Ill A we present and discuss exact results for the interaction energy Vr[7jj] of the 
dimerized Hubbard model on representative finite and infinite chains. These are obtained 



by solving Eq. (16) using accurate numerical methods. The dependence of the interaction 
energy on the alternating NN density-matrix elements 712 and 723 is analyzed. Scaling 
properties are identified within the domain of representability of jij. On the basis of these 
results we propose, in Sec. [Ill B| , a simple general approximation to ^^(712,723) which is 
useful for practical applications. A first test on the accuracy of this approximation is also 
provided by comparison with available exact solutions. 



A. Exact calculated VF[7ij] of the dimerized Hubbard model 

In the following we consider the dimerized ID Hubbard model which in the usual notation 
is given byi 

H =J2 t^j4aCja + U ^H^i] ■ (17) 
{ij)a i 

The NN hopping integrals tij take two alternating values: = ti2 = t + 5t for i odd and 
tj^j+i = = t ~ 6t for i even. The corresponding interaction-energy functional reads 



Vr[7ij] = min 



(18) 



where the minimization is performed with respect to all A^-particle states |\E'[7jj]) satisfying 
[7jj] I X^cr clo-Cjo- 1^1/ [7ij]) = 7ij for NN ij. For repulsive interactions M/^[7ij] represents the 
minimum average number of double occupations corresponding to a given degree of electron 
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delocalization, i.e., to a given 7jj. Eq. (|T^) then reduces to 

J2 A*, clcj^ 1^) + UJ2 ^i^f^n \^)=e\^) . (19) 

(ij) i 
a 

This eigenvalue problem can be solved numerically for finite systems with various boundary 
conditions. To this aim we expand |\l/[7jj]) in a complete set of basis states |$m) which have 
definite occupation numbers at all orbitals icr: njo-|$m) = ^i^l^m) with = or 1. 
The values of satisfy the usual conservation of the number of electrons A^e = + ^ej, 
and of the z component of the total spin Sz = {Ne^ — Nei)/2, where Ne^ = For 
not too large clusters, the state |^E'o[7jj]) corresponding to the minimum in Eq. (0) — the 
ground state of Eq. ( [T9| ) — can be determined by sparse-matrix diagonalization procedures 
such as the Lanczos iterative method.il For large chains, the properties of |\l/o[7ij]) can be 
calculated using the density-matrix renormalization-group (DMRG) methodil which allows 
reliable extrapolations to the infinite- length limit. Finally, in the absence of dimerization 
{6t = 0) translational symmetry implies that all NN 'jij are the same, and therefore one may 
set Xij = A for all NN ij. The lowest eigenvalue of Eq. (|l^) can then be determined from 
Lieb and Wu's exact solution of the ID Hubbard model following the work by Shiba.il 

In Fig. |l| the interaction energy W of dimerized Hubbard chains is shown in the form 
of constant-energy curves given by ^^(712, 723) = XE^f, where 712 = 7ii+i for i odd and 
723 = 7i,i+i for i even are the density-matrix elements or bond orders between NN's. -Ehf = 
U/A stands for the Hartree-Fock energy, and A is a constant (0 < A < 1). Results are 
presented for the A^^^ = 12 site ring and for the infinite ID chain which were obtained 
from Eq. (p!9| ) by using Lanczos-diagonalization and DMRG methods, respectively, i'i Only 
positive 7i2 and 723 are considered since this is the relevant domain when all the hopping 
integrals have the same sign. In bipartite lattices, like open chains or rings with even Na, 
the sign of the NN bond orders can be changed without altering W by changing the sign of 
the local orbitals at one of the sublattices. Thus, 1^(712,723) = W^(— 712, —723)- Moreover, 
^^(712,723) = ^^(723,712) as even an odd sites can be interchanged by a simple translation 
{Na even for rings). 

The domain of definition of W is restricted by the pure-state representability of jij. The 
axes 7i2 = and 723 = in Fig. |l| represent a collection of disconnected dimers or fully 
dimerized states, while 712 = 723 corresponds to non-dimerized states. In between, the 
degree of dimerization can be measured by the angle = arctan(7i2 /723). The degree of 



electron delocalization for each (f) is characterized by 7 = y7i2 + 723; which is bounded by 
7°°(0) < 7 < 7°(0) in order that 'jij remains pure-state representable. The density-matrix 
elements along the curve 7 = 7°(0) are the largest bond orders that can be achieved on a 
given lattice and for given Na and A^e [(7i2!723) = 7°(0)(cos0, sin0)]. They represent the 
maximum electron delocalization for each and yield the extremes of the kinetic energy 
Ek = J^tijlij, with different corresponding to different ^12/^23- Thus, for 7 = 7°(0) 
the density matrix can be represented by the ground state of the uncorrelated Hubbard 
model for some tu/hs iJJ = 0). In the absence of degeneracy the underlying electronic 
state |\E'o) is a single Slater determinant and W^(7i2 ) 723) = -^hf- Consequently, the upper 
bound for 7 coincides with the A = 1 curve in Fig. |l|. The correlation energy Ec = W — E^p 
vanishes as expected in the fully delocalized limit. For U = the minimization of the energy 
E = Ex as a function of 7jj can be stated in terms of the representability of 7jj alone. In 
this case the equilibrium condition yielding 'jfj is achieved at the borders of the domain 
of representability, more precisely, when the normal to the curve 7 = 7*^(0) is parallel to 

— * 

^Ek = {ti2, ^23)- 

Concerning the lower bound 7°°(0), one should first note that as 7 decreases, 7 < 'y^{4>), 
it is possible to construct correlated states |^E'[7ij]) having increasingly localized electrons. 
Charge fluctuations can then be reduced more efficiently for smaller 7, and therefore the 
Coulomb interaction energy decreases with decreasing 7 [see Eq. (|18|) and Fig. W reaches 
its minimum value Woo = Umax{0, — Na} in the strongly correlated limit where 7 = 
7°°(0). For half-band filling this corresponds to a fully localized state having 7°°(0) = 
and Woo = 0. However, note that for Nf, 7^ A^^^, W reaches W^o already for 7°°(0) > since 
partially delocalized states can be found having minimal Coulomb repulsion. This is the 
case for example in a fully-polarized ferromagnetic state. 

Fig. |l| also provides a qualitative picture of the functional dependence of W for dimerized 
chains. On the one side, for strongly dimerized states (0 ~ or ~ 7r/2) the constant- 
ly curves resemble circumference arcs, the gradient VW being approximately radial. This 
type of behavior is most clearly seen for weak or moderate correlations (A > 0.3), while 
in the localized regime (A < 0.1) it holds only for a very limited range of around = 
or = 7r/2. On the other side, for weakly to moderately dimerized states (0 ~ 7r/8- 
7r/4) the level curves can be regarded in first approximation as straight lines parallel to 
712 = —723- The very weak dependence of VW on implies that for (p c:^ tt/A the ground- 



10 



state values of 7^2 and 723, which result from the minimization of = Ex + W, are very 
sensitive to the hopping alternation 6t. In fact, significant variations of (p are necessary until 
VW = —VEk oc {\ + 5t/t^l — 5t/t) even ioi 5t/t 1. This is particularly notable for weak 
correlations since the W = E-^p curve is strictly linear for |0 — 7r/4| < 0.05. Therefore, a 
discontinuous change from 712/723 = 1 to 712/723 = 0.91 is found at f/ = and arbitrary 
small 6t. For U > 0, 712 and 723 are continuous functions of 6t, although the dependence 
on 5t remains very strong for small 6t, as can be inferred from the level curves in the figure. 
Comparing subfigures (a) and (b) one observes that the results for Na = 12 and Na = 00 are 
quite similar. The rather rapid convergence with chain length suggests that 11^(712,723) is 
not very sensitive to the details of the considered system, even if the minimization constraints 
in Eq. (p!8D apply only to NN bond orders. This is of interest for practical applications, as 
it will be discussed below. 

In Fig. ^ the interaction energy W is shown as a function of 7 for representative values 
of 0, including in particular the fuUy-dimerized (0 = 0) and non-dimerized (0 = 7r/4) cases. 
Despite the quantitative differences among the various 0, several qualitative properties are 
shared by all the curves: (i) As already discussed, the domain of representability of 7 is bound 
for each by the bond order 7°(0) in the uncorrelated limit. 7*^ decreases monotonously 
with increasing for < < 7r/4 showing that a compromise between 712 and 723 is made 
when the two bonds are active. This is an important contribution to the 0-dependence of 
W. (ii) In the delocalized limit, 1^(7°, 0) = Euf = U/A for all 0, since the electronic state 
yielding the largest 7 is a single Slater determinant. Moreover, one observes that dW/d'y 
diverges at 7 = 7° for all 0. This is a necessary condition in order that the ground-state 
density matrix satisfies 7^'^ < 7° for arbitrary small [/ > 0, as expected from perturbation 
theory, (iii) Starting from 7 = 7°, decreases with decreasing 7 reaching its lowest possible 
value ly = for 7 = {N^ = Na). The decrease of W with decreasing 7 means that the 
reduction of the Coulomb energy due to correlations is done at the expense of kinetic energy 
or electron delocalization. (iv) In the limit of small 7 one observes that ly oc 7^. Therefore, 
for t//t ^ 1, 7^'^ oct/U and Egs oc t'^/U, a well known result in the Heisenberg limit of the 
half-filled Hubbard model.0 
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B. Scaling approximation to VF[7ij] 

In order to compare the 7-dependence of W for different and to analyze its scaling 
behavior it is useful to bring the domains of representability for different to a common 
range by considering 14^(7,0) as a function of 7/7°(0), as displayed in Fig. ^. In this form 
the results for different appear as remarkably similar, showing that the largest part of the 
dependence of W on the ratio 712/723 comes from the domain of representability of 'jij given 
by its upper bound 7°(0) [7°°(0) = for half-band filling]. An analogous scaling behavior 
has been found in previous numerical studies of ^^(712) of non-dimerized Hubbard models, 
where 712 refers to the NN density-matrix element.0 In this case one observes that Vr(7i2) 
depends weakly on system size Na, band-filling n = N^/Na, and lattice structure, if W is 
measured in units of the Hartree-Fock energy E^p and if 712 is scaled within the relevant 
domain of representability [71^,712]- In the present context. Fig. ^ implies that the change 
in W associated to a given change in the degree of delocalization 7/7" (0) can be regarded 
as nearly independent of and system size. A good general approximation to W^(7, 0) can 
then be obtained by applying such a scaling to the functional dependence extracted from 
a simple reference system. An appropriate choice is provided by the fuUy-dimerized chain 
corresponding to = 0, which can be worked out analytically. In this case the system 
consists of a collection of dimers and the exact interaction energy reads 



Scaling the functional dependence of the dimer interaction energy to the 0-dependent domain 
of representability one obtains 



which we propose as approximation to W for dimerized systems. Notice that VFo(7; 0) 
preserves the previous general properties (i)-(iv) and that it is of course exact for = 
[7*^(0 = 0) = 1]. In practice, the system specific function 7°(0) can be easily obtained by 
integration of the single-particle spectrum. 

It is important to remark that the density matrices 'jij involved in the approximate 
functional Wq are pure-state iV-representable. Eq. (^) applies to the •jij obtained by scaling 
the off-diagonal elements of the density matrices 7°- that derive from uncorrelated states 




(20) 
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|\I/o) having Ne electrons on Na sites, and a uniform density distribution (\I/o| Z^o-^io-l^o) = 
Ne/Na = 1. In other terms, jij has the form jij = X'-ffj with < A < 1 for all i ^ j, 
and = 7° = 1 for all i. In order to show the pure-state representability of 7jj, we 
consider two normahzed A^-particles states |\E'a) and |\E';,) satisfying (\l'(j| X^o- c|o.Cjo-|^b) = 
for all ij. This condition is fulfilled, for example, by states a and b having different defined 
total spins S or Sz, or by superpositions of pure-S* or pure- 6*2 states sharing no common 
eigenvalues. The density-matrix represented by |^) = a;|\&a) +/5|\E'f,) with a'^+jS'^ = 1 is then 
given by 'jij = {"^l Z^o- cJo-Cjo-l^) = «^7j°- + P'^lijy where 7^" and 'y^j are the density matrices 
corresponding to |\E'a) and I^E'f,). Therefore, all the density matrices in the segment defined 
by 7j" and 'y'^j are pure-state A^-representable. The representability of a scaled uncorrelated 
7°- at half-band filling follows from the previous lemma by taking |\E'a) = I^E^o); which has 
5* = or 1/2, and |\E'b) equal to the fully localized state with one electron per site and 
maximal S = Ne/2, for which 7^^- = for all i ^ j, and 7fj = 1 for all i. Consequently, the 
Jij in the domain of definition of Wq and the ground-state density matrices jfj derived from 
it are all pure-state A^-representable. 

Fig. ^ compares Eq. (pT]) with the exact W{'y,(f)) for a 12-site Hubbard ring and for 
the infinite chain. One observes that the proposed approximation follows rather closely 
the exact results for all 7 and 0. The largest discrepancies are found for vanishing or 
moderate dimerization (e.g., (p = 37r/16 or = tt/4) and relatively large 7 (7 ~ 0.8). In 
all cases the quantitative differences remain small {\Wq — W\/U < 0.047 for = 7r/4 and 
\Wo — W\/U < 0.045 for = Stt/S) which is quite remarkable taking into account the 
simplicity of the approximation. In the following, Eq. (^If) is applied in the framework of 
LDFT to determine several properties of the dimerized ID Hubbard model. Comparison is 
made with exact results whenever possible in order to assess the performance of the method. 

IV. DIMERIZED HUBBARD CHAINS 

In Figs. I and |^ the ground-state energy Egg, kinetic energy E^, and Coulomb energy 
Ec of the ID Hubbard model are given as a function of the Coulomb repulsion strength U/t 
for different hopping alternations 6t. Accurate numerical results are also shown which were 
obtained using the Lanczos diagonalization methodil for A^^^ = 12 or the DMRG methodH 
for the infinite chain. In the case of the non-dimerized infinite chain, Lieb and Wu's exact 
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solutionE3 is taken as reference. The results for Na = 12 and Na = oo are qualitative very 
similar. Egg increases monotonically with U/t since dEgg/dU = {ni^fiii) > 0, vanishing 
in the limit of U/t = oo. For U/t < 4 this is essentially a consequence of the increase of 
Ec oc U, as Ex and jij remain very much like in the uncorrelated U = state. In contrast, 
for U/t > 4 the electrons become increasingly localized, and the increase of Eg^ results form 
the increase of Ek which approaches zero as decreases. At the same time Ec tends to 
zero as charge fluctuations are suppressed (see Figs. |]andH). 

Comparison between LDFT and the exact results shows a very good agreement. This 
concerns not only Egg but also the separate kinetic and Coulomb contributions indicating 
that electron localization and intra-atomic correlations are correctly described for all U/t. 
Moreover, this also shows that the results obtained for the ground-state energy are not the 
consequence of strong compensations of errors. Concerning the accuracy of Ek and Ec one 
generally observes that a somewhat higher precision is achieved for E^, which functional 
dependence is known exactly, as compared to Ec, which derives from an approximation 
to W [Eq. ([21|)]. For 6t/t > 0.1 the LDFT calculations are nearly indistinguishable from 
the exact ones (e.g., \Egs — Eg^\/t < 0.03 for 6t/t = 0.1). Even the largest quantitative 
discrepancies, found for the non-dimerized chain at intermediate U/t, are pretty small. For 
instance, for 6t = and U/t = 4 we obtain \Egs — Eg^\/t = 0.020 for the 12-site ring and 
\Egs — Eg^\/t = 0.044 for the infinite chain. Comparing Figs. ^ and ^ one observes that 
the performance of the method is sometimes higher for the 12-site ring than for the infinite 
chain. For example. Fig. ^ shows that Ek (Ec) is slightly overestimated (underestimated) 
for 6t = and U/{U + 4t) = 0.7-0.8, whereas for A^^a = 12 a much better agreement with 
the exact result is found (see Fig. |^). In any case it is important to recall that no artificial 
symmetry breaking is required to describe correlation-induced localization correctly, as it 
often occurs in other approaches (e.g., mean-field spin-density- wave state). Moreover, the 
present calculations remain simple and numerically not demanding, since the minimization 
of E[yij] is performed using analytical expressions for Ek and W [see Eqs. (H) and (pTI)]. One 
concludes that LDFT, combined with Eq. (^) as approximation to the interaction energy 
functional, provides an efficient and correct description of the ground-state properties of the 
ID Hubbard model in the complete range of interaction strength and dimerization. 
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The charge-excitation or band gap 

AE, = Egs{Ne + 1) + Egs{Ne - 1) - 2Egs{Ne) (22) 

is a property of considerable interest in strongly correlated systems which can can be related 
to the discontinuity in the derivative of the kinetic and correlation energies per site with 
respect to the electron density n. The determination of AE^. constitutes a much more serious 
challenge than the calculation of ground-state properties like Egs^ Ek, and Ec particularly 
in the framework of a density-functional formalism. Results for AEc of the ID Hubbard 
model are given in Figs. |^ and ^ as a function of the Coulomb repulsion strength U/t for 
different values of the hopping alternation 6t {n = 1). AEc vanishes for 6t = and U/t = 0, 
and increases with increasing U/t or 6t. Comparison between LDFT and Lanczos exact 
diagonalizations {Na = 12) or the Bethe-Ansatz solutionil [Na = oo and 6t = 0) shows 
fairly small quantitative discrepancies. In the most difficult non-dimerized case we find 
\AEc - AE^/l < O.m for Na = 12, and \AEc - AE^f\ < OMt for Na = oo. For small U/t 
and 6t = 0, AEc is somewhat underestimated for A^^^ = 12 and overestimated for Na = oo. 
The latter is mainly due to the fact that Eq. ( |2T| ) fails to reproduce the exponential decrease 
of AEc for U/t ^ {Na = oo and 5t = 0).@ As in previous properties the accuracy improves 
with increasing 6t. Fig. |^ shows that the LDFT results for non- vanishing dimerization and 
Na = 12 are very close to the exact ones {\AEc — AE^^\/t < 0.011 already for 6t/t = 0.1). 
Therefore, one expects that the predictions for Na = oo and 6t > should be reliable. 
Finally, one may note that in the limit of large U/t, the hopping alternation 6t has little 
effect on the charge gap. As the electrons tend to localize for f//t ^ oo, AEc ^ U + E^ 
where Eb = —^t is the energy of the bottom of the single-particle band. The present lattice- 
density-functional scheme describes correctly the crossover from a band insulator to a Mott 
insulator which occurs in dimerized chains as f//t is varied from the weak-interaction to the 
strong-interaction regime. 

V. SUMMARY AND OUTLOOK 

A novel density-functional approach to lattice-fermion models has been applied to the 
dimerized ID Hubbard Hamiltonian. In this framework the basic variable is the single- 
particle density matrix and the key unknown is the interaction-energy functional 
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VF[7ij]. In the present paper we have first investigated the functional dependence of W 
on the density-matrix elements 712 and 723 between nearest neighbors in dimerized chains 
{11,1+1 = 7i2 for i odd and 7i,i+i = 723 for i even). Rigorous results for VF(7i2 , 723) were 
derived from finite-ring Lanczos diagonalizations and from DMRG calculations for the in- 
finite chain. An analysis of these exact results shows that W can be appropriately scaled 
as a function 7/7^(0), where 7 = ^ 7^2 + 723 and 7°(0) is the largest representable 7 for a 
given = arctan(7i2/723). A simple general approximation to W was then proposed which 
takes advantage of this scaling behavior and which provides with a unified description of 
correlations from weak to strong coupling regimes. Finally, using this approximation, several 
ground-state properties and the charge-excitation gap of dimerized chains have been deter- 
mined successfully as a function of Coulomb repulsion strength U/t and hopping alternation 
St. 

The accuracy of the results encourages more or less straightforward apphcations of the 
present approach to related problems such as multi-leg ladders or the two-dimensional square 
lattice with first and second NN hoppings {t-t' Hubbard model). Moreover, the possibihty 
of generalizing the present scaling approximation to an arbitrary number of independent 
variables 7ij deserves to be investigated in detail, since it would open the way to applications 
in very low symmetry situations including metal clusters and disordered systems. 
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Tl2 

FIG. 1: Constant interaction-energy curves of the one-dimensional (ID) Hubbard model as given by 
(712, 723) = AE'hf, where £^hf = U/A is the Hartree-Fock energy and A a constant (0 < A < 1). 
The NN density-matrix elements are 7i,i+i = 712 for i odd, and 7i,i+i = 723 for i even. Results are 
given for (a) the Na = 12 site ring and (b) the infinite chain, both at half-band filUng (ATg = Ng)- 
A = 1 corresponds to uncorrelated states and defines the limit of representability of 7jj. Unless 
indicated the difference in A between contiguous curves is AA = 0.1. 
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FIG. 2: Interaction-energy W of the ID Hubbard model at half-band filling {Ne = Na) as a 
function of 7 = 7i2 + tIs for different values of = arctan(7i2 /723): (a) ring with Na = 12 sites 
(b) infinite chain. The density-matrix elements are 7i,i+i = 712 for i odd, and 7i,i+i = 723 for i 
even. 
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FIG. 3: Interaction-energy W of the ID Hubbard model as a function of 7/7'' for different (j) = 
arctan(7i2/723)- 7 = \J 7i2 + 723 7^(0) is the largest representable value of 7 for the given (j), 
which corresponds to the uncorrelated limit (0 < 7 < 7^^, see Fig. Results are shown for (a) the 
Na = 12 site ring and (b) the infinite chain, both at half-band filling. 
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FIG. 4: Comparison between the exact interaction-energy functional Wex of the Hubbard model 
and the approximation Wq given by Eq. (^). Results are given for the Na = 12 site ring (dashed) 



and the infinite chain (solid) as a function of 7 = \/7i2 + 723 different (p = arctan(7i2 /723)- 
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U/(U+4t) 

FIG. 5: Ground-state energy Egg = Ek + Ec, kinetic energy Ek, and Coulomb energy Ec of 
dimerized Hubbard rings with hopping integrals tij = t{l it 5t), Coulomb interaction U, Na = 12 
sites, and A'^e = Na electrons. The symbols are obtained from exact Lanczos diagonalizationsii 
and the solid curves correspond to the present lattice density-functional theory. 
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U/(U+4t) 

FIG. 6: Ground-state energy Egg = Ek + Ec, kinetic energy Ek, and Coulomb energy Ec of 
dimerized infinite Hubbard chains with hopping integrals tij = t{l±6t). The symbols are obtained 
using the density-matrix renormalization group method@ and the solid curves correspond to the 
present lattice density- functional theory (see Fig. ||). 
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FIG. 7: Charge-excitation gap AEc of dimerized Hubbard rings with Na = 12 sites, hopping 
integrals tij = t{l it 5t), and band filling n = N^/Na = 1. The symbols refer to exact numerical 
diagonahzation Ji and the curves to the present lattice density-functional approach. 
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FIG. 8: Charge-excitation gap AEc of dimerized ID Hubbard chains with hopping integrals tij = 
t{l lb 5t) and band filhng n = N^/Na = 1. The crosses refer to exact Bethe-Anstaz results for 
6t = (see also Fig. 0) . 
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